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ABSTRACT 

In recent years it has been shown that microlensing is a powerful tool for ex- 
amining the atmospheres of stars in the Galactic Bulge and Magellanic Clouds. The 
high gradient of magnification across the source during both small impact parameter 
events and caustic crossings offers a unique opportunity for determining the surface 
brightness profile of the source. Furthermore, models indicate that these events can 
also provide an appreciable polarisation signal - arising from differential magnifi- 
cation across the otherwise symmetric source. Earlier work has addressed the signal 
from a scattering photosphere for both point mass lenses and caustic crossings. In a 
previous paper, polarimetric variations from point lensing of a circumstellar envelope 
were considered, as would be suitable for an extended envelope around a red giant. 
In this work we examine the polarisation in the context of caustic crossing events, the 
scenario which represents the most easily accessible situation for actually observing 
a polarisation signal in Galactic microlensing. Furthermore we present an analysis 
of the effectiveness of using the polarimetric data to determine the envelope proper- 
ties, illustrating the potential of employing polarimetry in addition to photometry and 
spectroscopy with microlensing follow-up campaigns. 
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1 INTRODUCTION 

Much attention has been paid to the situation in which the 
analytic case of the magnification of a point source by a 
point lens breaks down. It has been noted that for small lens- 
source separations, the finite size of the source star needs to 
be considered (Nemiroff & Wickramasinghe 1994; Witt & 
Mao 1994; Witt 1995; Peng 1997). Such events not only con- 
strain the lens properties by breaking the degeneracy in the 
event parameters (Gould 1994) but also provide valuable stel- 
lar atmosphere information, such as limb darkening (Valls- 
Gabaud 1998; Hendry et al. 1998), polarisation (Simmons et 
al. 1995a,b; Newsam et al. 1998), motions in circumstellar en- 
velopes (Ignace & Hendry 1999) and the presence of starspots 
(Heyrovsky & Sasselov 2000). The opportunity for study- 
ing stellar atmospheres through microlensing as described in 
Sackett (2001) and Gould (2001) has clear advantages over 
other methods such as eclipsing binaries, because the source 



and means of studying the source are not coupled and the flux 
from the source is magnified rather than diminished. 

Binary lenses produce another deviation from a standard 
microlensing lightcurve, because caustics are produced (e.g., 
Schneider & Weiss 1986). It is the high gradient of magnifi- 
cation across the caustic that allows one to infer information 
about the source intensity profile, even though the source im- 
ages are not individually resolvable with current instruments. 
Despite binary lens events only accounting for about 5% of 
microlensing events (e.g., Alcock et al. 2000), it is more likely 
that finite source effects will be relevant for a binary lens rather 
than a point lens. The reason is that finite source effects are 
mainly discernible only when the lens transits the source it- 
self (Gould 1994). Since the angular Einstein radius (see 
Eq.|3j is usually much larger than the angular source size, such 
transits tend to be rare. On the other hand, the caustic struc- 
tures giving rise to high magnifications in binary lens events 
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are spatially extended (of order 0b in scale), and so caustic 
crossings by the more distant source are relatively common in 
events associated with binary lenses. So it is more expedient to 
study the resolution of stellar atmospheres due to binary lenses 
than point lenses. The structure of the caustics produced and 
the resulting magnification and lightcurves from binary mi- 
crolenses has been discussed extensively in the literature (Mao 
& Paczyhski 1991; Mao & Di Stefano 1995; Di Stefano & 
Perna 1997; Dominik 1998; Gaudi & Gould 1999; Dominik 
2004b). A number of papers have reported on source prop- 
erties (such as limb darkening) derived from the light curves 
of binary lens events (Albrow et al. 1999; Afonso et al. 2000; 
Albrow et al. 2000; Albrow et al. 2001a; Albrow et al. 2001b; 
Abe et al. 2003; Cassan et al. 2004). 

Considerable theoretical work on the polarisation signa- 
tures from microlensing has also been carried out. Schneider 
& Weiss (1986) were the first to discuss the use of caustic 
crossings for inferring source intensity profiles. Schneider & 
Wagoner (1987) calculated the polarisation from the lensing 
of thick scattering photospheres of supemovae during such 
crossings. Simmons et al. (1995a,b) examined the polarisation 
from an electron scattering atmosphere being microlensed by 
a point mass, thus allowing the limb polarisation to be mea- 
sured. Newsam et al. (1998) used this analysis to demonstrate 
that even relatively 'poor' polarimetric data can considerably 
improve the determination of stellar radii. In turn, Agol (1996) 
modelled the polarisation from an electron scattering stellar 
atmosphere by a binary lens. More recently, Simmons et al. 
(2002) examined the polarisation from an extended envelope 
for a point mass lens. The work of this paper applies the same 
atmosphere model as in Simmons et al. (2002) to the case of 
caustic crossings - arguably the 'best case scenario' for prob- 
ing the extended envelope. 

The structure of this paper is as follows. In Section 2 we 
discuss the microlensing of extended sources. In Section 3 we 
describe the polarisation intensity maps used in the calculation 
of the microlensing lightcurves. In Section 4 we provide repre- 
sentative lightcurves for both the polarisation and flux signals 
for a range of source parameters. We also discuss the obser- 
vational implications of our results. The duration of caustic 
crossing events typically occurs over one night, so we believe 
that it is important to explore different observing strategies to 
find the most suitable way of observing these events while ef- 
fectively recovering the source and envelope properties rela- 
tive to the Einstein radius. Concluding remarks are presented 
in Section 5. 



2 MICROLENSING OF EXTENDED SOURCES 

The magnification of a point source by a point lens is given by 



(1) 



where u is the angular separation between the lens and source 
as normalized by the angular Einstein radius 6'e, which is 
given by 



where Ml, Ds, and Dl are the lens mass, the distance to 
source and the distance to lens respectively. The lensing is a 
transient event with u = u(t) a function of time given by 



u{t) = 



tl 



(3) 



where to is the time of maximum magnification at ito the im- 
pact parameter. The parameter tE is the crossing time of 6^: 



(4) 



Se = 
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where /irci is the relative proper motion between the lens and 
the source. 

The point source approximation is valid when 9s <^ Oe, 
for 6s the angular radius of the source (which we shall assume 
is circularly symmetric). When this no longer holds - such as 
when the lens transits the source - the magnification is de- 
termined by an intensity-weighted integral over the projected 
surface of the source. The resultant is, 

Anct(t) — „ g ^ (3) 

J^^ i{e,a)ededa 

where 6 is an angular radius measured from the source ori- 
gin, and a is an azimuthal angle about that origin. The angle 
9l is the angular separation between the point lens and the 
source centre. The integral means that different emitting el- 
ements of the projected source make weighted contributions 
to the microlensing light curve according to their relative pro- 
jected proximity to the lensing mass. Since this proximity is a 
function of time, an analysis of the event light curve allows for 
the possibility of determining the surface brightness profile of 
the source. 

The magnification due to a binary lens system depends on 
additional parameters relating to the binary system: the mass 
ratio, the separation of the lenses, and the angle defining the 
trajectory of the source relative to the binary lens orientation. 
The magnification due to a binary lens has no simple analytic 
dependence on these parameters. However, near the fold caus- 
tic, the magnification can be approximated by 

Acau=^0 + -^, (6) 

where d — d(t) is the angular distance (normalized to 6e) 
from the caustic to a source element for all such elements in- 
terior to the caustic, Aq is the constant magnification of the 
three non-caustic images, and &o is a scale factor related to 
specifics of the lens geometry and the ratio of 6*3 to 6'e (e.g., 
the discussion of Castro et al. 2001). The literature of reports 
on caustic crossings by stars indicate typical observed peak 
magnifications of around 20-30 are common (e.g., Abe et al. 
2003; Afonso et al. 2000; Albrow et al. 1999; Albrow et al. 
2000; Albrow etal. 2001a; Albrow et al. 2001b). The sources 
are generally red giant stars, with radii of around R, ~ lOJ?©. 
With es/fe ~ 10^ values of 60 are around 10, and values 
of Ao are of order 3-10, as evidenced by the 'troughs' be- 
tween successive caustic crossings (e.g., Alcock et al. 2000). 
The particular form of the approximation in Eq. |6| assumes 
that the caustic is a straight line (i.e., the source must be small 
compared to the curvature of the caustic) and that the caustic 
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crossing point is not in the vicinity of a cusp, wliere the magni- 
fication function takes a different form. Exterior to the caustic 
(yet near it so that the curvature of the caustic can still be ig- 
nored), the magnification will be constant with Acau ~ Aq. 
Note that the time-dependent variation of the polarised flux 
results entirely from the second term of Eq.|6|for intrinsically 
symmetric and unresolved sources. 



how the scattering opacity evolves in the flow. We additionally 
assume the envelope is optically thin for this application, but 
will explore optical depth effects of the envelope in a subse- 
quent paper. 

We adopt the Stokes vector notation I = (/, Q, U, V)''" 
for the intensities, and F = (Fj, Fq , Fu , Fv ) for the fluxes. 
Following the notation of Simmons et al. (2002), the Stokes 
intensities for direct and scattered light will be given by 



3 POLARIZATION MODEL 

We review now the intensity profile that will be used with Eq. 
|5|to calculate the microlensing lightcurves. As stated above 
the formalism follows that of Simmons et al. (2002). As em- 
phasized in that work the model employed is well-suited to 
evolved cool stars. This class of stars exhibits stellar winds 
that are significantly stronger than the Sun's, with mass-loss 
rates ranging from 10"^" A/q yr~^ for typical red giants up to 
1O~'''M0 yr~^ for red supergiant and asymptotic giant branch 
stars (e.g., Earners & Cassinelli 1999). The more extreme stel- 
lar winds are clearly dust-driven (e.g., Netzer & Elitzur 1993; 
Habing, Tignon, & Tielens 1994). Stars with milder winds are 
less understood, possibly driven by a combination of a Parker- 
style wind augmented by molecular and dust opacities (Jor- 
gensen & Johnson 1992). 

Assuming the flows are spherically symmetric, the run of 
the bulk gas density p in the wind with radius r will be given 
by 



P = 



A'/ 



Anr^ v{r) ' 



(7) 



where Af is the mass-loss rate, and v{r) is the radial veloc- 
ity of the flow, beginning subsonically and asymptoting to a 
terminal speed Voo at large radius. 

The scattering opacity responsible for producing polar- 
isation could be molecular Rayleigh scattering or dust scat- 
tering. Consequently, the density profile of the polarigenic 
species can differ from Eq. by virtue of how molecules 
or dust are produced and destroyed as a function of radius. 
Particularly in the case of dust, the location of the conden- 
sation radius will be important. This latter point concerning 
dust is interesting, since the stars under consideration have 
photospheric effective temperatures of 3000-4000 K, whereas 
the dust condensation temperature is typically around 1500 K 
(Gail & Sedlmayr 1986). The implication is that dust-driven 
winds can have central cavities (not vacuums, but rather inte- 
rior zones for which there is no significant scattering opacity). 
Models for dust-driven winds naturally predict the radial ex- 
tent of these cavities, yet the condensation radius is largely un- 
constrained by observations (Bloemhof & Danen 1995). Mi- 
crolensing can provide relevant observational constraints on 
the location of the condensation radius. For this pilot study 
of the polarimetric signals from caustic crossing events, we 
choose to parametrize the scattering number density by a sim- 
ple power-law with 



n{r) = no (^) 
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where _Rh is the radius for any 'hole' of scattering opacity, and 
the total optical depth, Tsc is defined as 



(10) 



with a the scattering cross section, and where the scatterers 
are taken to exist only for r 7?h. 

For simplicity we shall assume a uniform surface bright- 
ness profile for the star, such that 



Mp) 



(/., 0,0,0) forp<i?* 
for p > R, 



(11) 



where /, is the intensity of the star. Although we ignore limb 
darkening, its main effect would be to enhance the polari- 
sation from lines-of-sight that are close to the star (bearing 
in mind that we are currently ignoring scattering polarisation 
from the stellar photosphere, an assumption that we shall jus- 
tify a posteriori). The reason is that at a fixed distance from 
the star, limb darkening causes the stellar radiation field im- 
pinging upon a given point at this distance to be more radially 
directed than is the case for a uniformly bright star (Cassinelli, 
Nordsieck, & Murison 1987). 

Also appearing in Eq. llOl is the stellar occultation factor, 
go, that is given by 



, , [1/2 forp < _R. 
^«(^) = ^ 1 for p>i?. 



(12) 



The occultation factor simply accounts for the fact that radi- 
ation scattered on the far side of the star will not reach the 
observer. The integral factors Gi and Gp are 



Gi = 

and 

Gp = 



1 - s 



3 y s'^q 



(13) 
1 Ids 



instead of fretting about the details of the wind acceleration, or 



1 — s 



with s = {p/r)^, q = (R-t/p)"^, Smax = (p/rmin)^, for 



(14) 
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rmin = max(i?h,p), 



(15) 



4.1 Test Case 



where z is the hne-of-sight coordinate and — + . 

The fluxes during the microlensing event are then com- 
puted from the integral expressions: 



/• oo /'27r 

Fi = / I{p,a)Ac^u{d)pdpda, (16) 

Jo Jo 

poo p2n 

Fq = / I{p,a) cos 2a Aciuu{d)p dp da, (17) 

Jo Jo 

Fij = / I{p,a) sin2a Acau(rf)pdpda, (18) 

Jo ^0 

Fv = 0. (19) 

As is usual, the observed fractional polarisation can then be 
calculated as 



P = 



^F^ + F-+F, 



Fi 



(20) 



where the total flux Fi is a sum of the direct stellar flux F, and 
the scattered intensity Fi^^ . The polarisation position angle is 
defined as 

7A = itan-^||. (21) 



4 RESULTS 

Having described the lensing approximation for a straight fold 
caustic, and the underlying source model, we have conducted a 
parameter study for microlensing light curves associated with 
caustic crossing events as various source and lens properties 
are varied. In displaying results, our goal is to highlight those 
features that pertain to elucidating the properties of the source; 
therefore, instead of light curves as a function of time, we 
choose to plot flux observables as a function of caustic position 
relative to the star. To do so, we define the coordinate aliens as 
the normal projected distance between the straight line caustic 
and star centre in the source plane. Then yiens is the coordinate 
along the caustic direction. Both coordinates are normalized to 
the stellar radius R, . The case aliens < is when the star cen- 
tre lies inside the caustic; the case xicns > is when the star 
centre lies outside the caustic; and a;ions = is the moment of 
transit for the star centre. 

Conversion to a time coordinate is achieved with t = 

(liens R*)/ {Ds Mrcl.x), whcrC ^roLx = Mrol COS 7, with /^rcl 

the magnitude of the relative proper motion between the 
source and caustic, and 7 is the trajectory orientation of the 
source relative to the caustic. So 7 = 0° means the source 
is traveling in the +x direction in the frame of the caustic, 
whereas 7 = ±90° means the source is moving in the ±y di- 
rection, respectively. The orientation of the trajectory has no 
bearing on the shape of the lensing light curves in time, except 
as a 'stretch' factor for scaling purposes; for example, 7 does 
not affect the value of the peak polarisation achieved during 
the event. 



Figure Q shows model results for a test scenario. The upper 
panel shows the flux of the source, relative to the case of being 
unlensed, as a function of the caustic location. The lower panel 
is for the polarised flux, also normalized to the unlensed flux 
of the source Fo, and multiplied by 100 to simulate a kind of 
percent polarised flux (the reason for this rather odd choice 
will become apparent in a moment). The envelope parameters 
are fixed with an envelope optical depth of Tsc =0.1, density 
parameter /3 = 2, and no cavity (i.e., Rh — R*). The lens 
parameter bo is fixed at a value of 10, but Ao is allowed to 
vary from 1 to 15, with larger values giving stronger peak flux 
magnifications. 

We have argued that Ao should have no bearing on the 
polarised emission, and indeed that is seen to be the case in 
FigureQ All of the variation of the polarized flux comes from 
the second term in Eq.|6|with scaling 60. In the standard rep- 
resentation of polarisation with p = Fq / Fi . The fractional or 
percent polarisation can be affected by the value Ao, because 
the total flux variation as the microlensing event evolves does 
depend on Ao- A value of Ao = 1 has been adopted for the 
rest of the model calculations. 

It is useful at this point to introduce a schematic figure 
that demarcates regions contributing to the polarised light as 
the caustic crossing evolves. Figure|2|shows four panels, with 
the source moving from inside the caustic to outside in the se- 
quence (A) to (D). The vertical dotted line is the caustic, hence 
the left region is interior to the caustic and the right region is 
outside it. The star is the cross-hatched region. The double- 
headed arrows show the sense of orientation that would result 
for the emergent polarised flux from right angle scattering in 
the plane of the sky if the source were resolved. The dashed di- 
agonal lines are where Q = 0. Thus the figure is useful in map- 
ping how the different zones of polarised flux will contribute to 
the total polarised emission from the unresolved lensed source 
as a function of its location relative to the caustic line. 

The variation of the model polarised flux shows inter- 
esting sign changes. Here we are assuming an observational 
scenario in which a second exiting caustic crossing has been 
predicted from an earlier interior crossing event. Thus in Fig- 
ure|2|at time (A), the source is approaching the caustic to exit. 
The quadrant closest to the fold is dominated by scattered light 
with Q > 0, and so initially the polarisation is oriented paral- 
lel to the fold caustic in the sky. As the event progresses, the 
situation of case (B) is reached. Only scattered light leftward 
of the caustic will be subject to a magnification gradient, thus 
breaking the symmetry and leading to a net observed polarisa- 
tion. By time (B), the light being most strongly magnified is 
for Q < 0, and so the polarisation position angle has rotated 
90° by the time the caustic is first tangent to the photosphere. 
By time (C), a minimum in the polarised flux has passed. Now 
the polarisation changes rapidly, so that by the time the caus- 
tic is tangent to the far side of the photosphere, the polarised 
flux has changed sign again, becoming positive. Note in terms 
of the polarised light Fg/Fi, the percent polarisation will rise 
significantly, because the magnification of the stellar photo- 
sphere is minimized once the photosphere has completely ex- 
ited the caustic. 
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Figure 1. Plots of the variation in total flux (upper panel) and the normalized polarised flux (lower panel) for a scattering envelope during a fold 
caustic crossing event. The lower axis displays the location of the caustic in the source plane (normalized to ft* ) with respect to the centre of the 
star, with negative values for when the star lies interior to the caustic, and positive values for when it lies exterior. Values of /3 = 2, = ij, (no 
cavity), and an envelope optical depth of Tsc = 0.1 are adopted. A lens parameter of 6o = 10 is fixed, but Aq takes values of 1, 3, 5, 10, and 15 
(in order of increasingly stronger peak fluxes). As described in the text, variations in Aq change the total flux light curves, but not the polarisation, 
because Aq is a constant of the integration on either side of the fold caustic. It is bo and the properties of the envelope that set the shape and 
amplitude of the polarimetric variations. 



4.2 Variable Envelope Optical Depth 

Figure |3| shows the response of the flux and polarisation light 
curves to different values of envelope optical depth. The lower 
panel now plots the standard form of polarisation, p = _Fq / Fi , 
and hereafter Values of r^c = 0.001,0.01,0.1,0.3, and 1.0 
are used. A value of Aq = 1 and 6o = 10 are used with 
a density distribution described by /3 = 2. These values for 
A(i,ho, and /3 will be standard in our calculations unless noted 
otherwise. The envelope has little effect on the flux light curve 
except when t^c > 0.1. On the other hand, the amplitude of 
the polarisation light curve is approximately linear in t^c- 

The envelope is reasonably thin, and so the flux magnifi- 
cation is dominated by the photospheric emission. Agol (1996) 
has investigated the polarimetric variations from scattering po- 
larisation in stellar atmospheres. We have purposely ignored 
photospheric contributions to the polarised emission, because 



circumstellar envelopes are more efficient at producing po- 
larised emission (albeit, this is a function of optical depth), 
and because we wish to investigate the effects of a circumstel- 
lar envelope for the light curves. Although Ao and ho are not 
exactly the same, the case of the solid line in the upper panel 
of Figure|3|is roughly comparable to the r = 0.01 case shown 
in the upper panel of Agol's Figure 1. 

Overall, the polarisation curves are generally similar to 
those of Agol (1996); however, there are notable quantitative 
and qualitative differences. First quantitatively, the peak po- 
larisation achieved by photosphere crossings were rarely in 
excess of 1%, and in some cases only a few tenths of a per- 
cent, whereas thin scattering envelopes can achieve values in 
excess of 5% when Tsc is large enough (tsc > 0.5). 

There are qualitative differences as well. The underlying 
source models are drastically different. For example, a stellar 
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Figure 2. Schematic of the geometi'y for the transit of the star and circumstellar envelope across a fold caustic. Shown are four times (A) - (D). 
The star is cross-hatched. The vertical dotted line is the caustic, with left being interior and right being exterior. The double-headed arrows are 
representative polarisation vectors at those points (vertical being Q > and horizontal being Q < 0). Finally, the diagonal dashed lines are the 
locus of points where the emergent Q-intensity would be zero; these then are 'null' lines where Q switches sign, and the polarisation position angle 
rotates 90° between adjacent zones. 



photosphere has an intrinsic polarisation profile (as a function 
of p) that is maximum at the stellar limb and decreases to zero 
at the centre of the star. In our case the photospheric polari- 
sation is ignored because it is small compared to the envelope 
polarisation. The extended envelope has a polarisation profile 
that is zero at large distance from the star, and increases toward 
the stellar limb. The polarisation peaks outside the stellar limb, 
and decreases to zero again at the stellar centre. So the source 
models are quite different, but at the same time the variable 
polarisation from microlensing for the two cases are some- 
what similar. Basically, the photospheric polarisation has a dis- 
continuous jump in moving from off the star across the stellar 
limb, whereas the circumstellar envelope has a more gradual 
peak off the stellar limb. Microlensing involves a weighted sur- 
face integral, thus 'smoothing' over these detailed differences, 
leading to somewhat similar light curves (Dominik 2004a). 
Still, the fact that the peak polarisation appears at the stellar 
limb for a photosphere in a discontinuous way explains why 



there is a peak (sometimes cuspy) at a::icns = —R* in Agol's 
models and not in ours. (This is = — r in Agol's notation; 
see his Fig. 1.) 

4.3 The Influence of a Cavity 

Figure|4|shows how a cavity of scatterers at the inner envelope 
impacts the polarisation light curves. The different curves are 
for different 'hole' radii of i?h = 1.0, 1.5, 3.0, 5.0 and 8.0i?.. 
Clearly, as the extent of the cavity increases, the polarimetric 
variation occurs over a longer time-scale, with the peak polar- 
isation shifting toward larger values of positive a:;icns, and the 
negative 'trough' growing in extent toward negative siens as a 
precursor to the transit of the star. 

In each case the envelope optical depth is maintained at 
Tsc = 0.1. First, this is rather optically thin, as evidenced by 
the total flux light curve (upper panel) which does not vary 
much between the different cases, and is dominated largely by 
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Figure 3. Plots like Fig.0 but now with fixed lens parameters and variable envelope optical depths, and with the lower panel as percent polarisation 
(normalized by the total intensity flux Fi). In this case values of /9 = 2, R]^ = _R* (no cavity), Aq = 1, bo = 10 are assumed. The different curves 
are for different envelope scattering optical depths, with Tsc = 0.001, 0.01, 0.1, 0.3, and 1.0 (in order of increasingly stronger peak polarisations). 
Note the change in sign of the polarisation as the stellar photosphere transits the caustic, and the strong peak value that results immediately after the 
photosphere has completely exited the caustic. 



the photosphere. Second, the scale of the polarisation is to ze- 
roth order determined by the value of Tsc, which explains why 
all of the curves have similar peak polarisation values. Third, 
holding Tsc constant implies conserving the total number of 
scatterers; hence although the different cases shown have the 
same density distribution at /3 = 2, these require different den- 
sity scales no, so as to maintain fixed Tsc- Using Eq. <10t . the 
density scale for fixed envelope optical depth is given by 



no = (/? - 1) 



(22) 



So our approach for inserting a cavity is not to delete scatter- 
ers, but to redistribute them outward. The goal of considering 
cavities is to illustrate how microlensing can neatly trace the 
cavity extent through the polarisation light curve, which is rel- 
evant for the case of red giants that can form dust in their winds 
at a condensation radius that is offset from the stellar photo- 
sphere. 

Figure |5| shows the same curves of Figure |4| but with 



xicns normalized to i?h instead of R, , which nicely shows how 
the polarimetric variations are set by the crossing of the cavity. 
Clearly, the peak polarisation is affected by the cavity (gener- 
ally smaller), and substructure is seen around the passage of 
the photosphere from inside the caustic to outside. The peak 
polarisation consistently occurs just before the cavity transits 
entirely out of the caustic, and so becomes an excellent tracer 
of the cavity extent relative to the stellar radius, which can be 
determined from the total flux variations. 



In other words the total flux variations show a peak at a 
time when the stellar limb just begins to transit the caustic. 
As the event proceeds, the total flux shows a precipitous drop, 
and then goes flat. The extended envelope can produce a tail of 
enhanced brightness after the star has transited out of the caus- 
tic, but the drop is dominated by the star. So that time-scale is 
2tf = 29,/firci, where 9, is the angular radius of the star. 
On the other hand, from this time until the peak polarisation 



8 R. Ignace, J. E. Bjorkman, and H. M. Bryce 




Figure 4. Illustration of how a cavity of scattering opacity impacts the polarisation variation. All of the curves are for Tsc =0.1. Each is distin- 
guished by the extent of the cavity, with R^^ = 1.0, 1.5, 3.0, 5.0, and 8.0-R», with larger cavities yielding peak polarisations at larger values of 
aliens - Note that the total flux curve (upper panel) is little influenced by the cavity extent because the scattering envelope is optically thin. Although 
basically similar, there are some notable qualitative differences in the variation of polarisation when a cavity is present (see Fig.lsl. 



is achieved will require a time = dh/l^rci- Consequently 
Rh/R* will be given by tt^/t^. 



4.4 The Impact of the Density Distribution 

Models for the flux and polarisation variations have been gen- 
erated for different envelope density distributions, and the 
results are displayed in Figure |6| The curves are for /3 = 
1.5, 2.0, 2.5, 3.0 and 4.0. A cavity with 7?h = 3.0R, has been 
adopted. 

The different cases are all for a fixed value of Tsc ~ 0.1. 
As noted before, this is in essence achieved via redistribution 
of scatterers. In this case there is a fixed hole. Changing /3 
makes the density distribution more or less steep. As /3 is made 
to increase, keeping Tsc fixed results in increased values of no, 
and so the peak polarisation that is dominated by the number 
density of scatterers at the limb of the cavity increases as well. 

So (3 does not necessarily lead to larger polarisations; 
here it is an artefact of maintaining a constant value of Tsc. 



Fundamentally what /3 does is to alter the slope of the polari- 
sation curves after the cavity has completely passed out of the 
caustic, making the slopes steeper with increasing /3. In fact, 
in the limit that the star can be treated as a point source, the 
asymptotic slope of the polarisation for relatively large values 
of aliens can be derived analytically. The derivation is found in 
the Appendix; here just the result is quoted. Asymptotically, 
the polarised flux (not the percent polarisation) will be given 
by 



i^Q OC 



' Rh 



(2/3-l)/2 



^-(2/3-l)/2^ 



(23) 



So larger values of /3 lead to steeper declines in the polarised 
flux as the microlensing event progresses. Two points should 
be mentioned. First, although real sources may not follow a 
power-law distribution for the density of scatterers in portions 
of their extended envelopes, a value of /3 = 2 is reasonable 
to expect at large scale for a spherical wind flow, for which 
case Fq oc t^^^^. Second, the preceding equation is only 
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Figure 5. The results shown in Fig.|4] except now the lower axis for the position of the fold caustic is normalized to the size of the cavity /J^. 
From Fig.|3 the total flux variation was set by the passage of the photosphere across the caustic. The polarimetric variations on the other hand 
are determined the size of the cavity. Although the positive peak polarisation has variable width, the negative 'trough' has approximately constant 
width. The solid line is for no cavity (i.e., iJh = R*)', notably the presence of a cavity changes the qualitative shape of the trough, with some 
recovery towai'd net zero polarisation followed by a sharp drop toward more negative values. 



valid both for when the star can be treated as a nearly point 
source of illumination with respect to the scattering envelope 
and when the caustic can be approximated as a straight line. 
Even if the asymptotic trend of Eq. <23> is not achieved in a 
real event, it does provide useful insight and limiting behavior 
for the modeling effort. 

4.5 Variation of the Magnification Gradient 

Figiire|7|shows model results as the value of bo is varied, with 
values of 3, 5, 8, 12, and 17 (with Aq = 1 fixed). With bo rel- 
atively large compared with Ao, the polarisation varies little 
for aliens < 1, prior to when the star has passed out of caustic. 
The reason is that the polarisation is a ratio of the polarised 
light to the total flux, and both scale as bo- Clearly, the to- 
tal flux shown in the upper panel is strongly dependent on 6o, 
increasing essentially linearly with bo - Similarly, the polarisa- 
tion after the star has passed out of the caustic is affected by 



bo, because now Fi is a constant that depends primarily on the 
photospheric flux multiplied by Ao, whereas Fq still depends 
on the value of bo. 



5 DfSCUSSION 

This study has demonstrated that polarimetric observations of 
caustic crossing events could be used to probe circumstellar 
envelopes. Such data could constrain the number density of the 
scatterers within the envelope, detect the presence and trace 
the extent of a central cavity around the source photosphere, 
and provide information about the density distribution of scat- 
terers. 

To summarize, a comparison of the time-scale for the to- 
tal flux variations against that of the polarimetric variations 
yields the extent of the cavity relative to the stellar radius. 
In the case of dust producing cool star winds, that informa- 
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Figure 6. Here only the value of /3 is allowed to vary. The flvix Ught cvirve hardly changes, whereas the polarimetric profiles, although similar 
qualitatively, are seen to vary significantly in amplitude. The values of /? are 1.5,2.0,2.5, 3.0, and 4.0, in order of stronger positive peak polarisations, 
and the envelope sports a cavity with ijji = 3.0-R,. 



tion is sufficient to test models that predict the location of the 
dust condensation radius. The late time evolution of the po- 
larised flux is set by the density distribution (our /3-value for 
this work). Lens parameters Ao and feo can be derived from 
fitting the variation of the photospheric flux during the caustic 
crossing, in which case model fits to the polarisation level will 
give the value of t^c for the circumstellar envelope. 

To connect our model results with applications to ob- 
served events, we must relate the envelope properties for the 
model to those of real sources. Some help toward this end is 
provided by Netzer & Elitzur (1993), who describe model re- 
sults for dust-driven winds. In their Eq. (10), they provide an 
expression relating stellar parameters M, Voo, and L» to the 
flux mean optical depth of the envelope tf: 

^ ^ _Wio^ 

2xlO-5M0yr-i Voo/lOkms"^' ^ ' 

here shown in slightly modified form from their paper. The 
value of TP will not equal the value of Tsc that we use to char- 
acterize our models; however, the flux mean opacity does give 



an overall scale related to the optical depth of the envelope. 

Although relating Tsc to rp will depend on the particular opac- 
ities involved, one might generally expect that the two will 
scale together. The models of Netzer & Elitzur show that the 
minimtmi mass-loss for dust driving to be dominant is around 
10"'' Mq yr"^ At this value for a star with L« = lO^'i© 
and Vac = 10 km s^^, the flux mean optical depth will be 
0.02 (of course, the optical depth at a wavelength of interest 
can be higher or smaller). 

One of the challenges in detecting polarized signals in 
real events, such as this hypothetical Bulge star, is that the 
crossing of the caustic by the photosphere will typically take 
only a few hours (e.g., Alcock et al. 2000). We can estimate 
the detectability of polarizations predicted by our models. Us- 
ing a Kurucz model for a cool subgiant of g = 3.5 and 
T = 4500 K (parameters similar to OGLE-1999-BUL-23 
from Albrow et al. 2001), the I-band flux at a distance of 8 
kpc is estimated to be 3 x 10~^^ erg s~^ cm~^, with a cor- 
responding magnitude of about mi « 18. Of course, during 
the lensing event, the source brightens, and magnifications by 
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Figure 7. These figures are for fixed envelope properties, but for different values of the lens parameter bo = 3, 5, 8, 12, and 17 and Aq fixed at 
unity. Both the strength of the peak flux and peak polarisation vary almost linearly with the value of 6o- Variations in Aq would have no effect on 
the polarisation, but would influence the total flux curves. 



an order of magnitude are typically achieved, at which point 
mj ~ 15.5. The time t required to achieve a given signal-to- 
noise ratio S/N for a telescope of diameter D in centimeters, 
allowing for Poisson noise only, will be 
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D2 



(25) 



Our models indicate that peak polarizations of about 1% will 
be achievable during caustic crossings. A 5a detection at 
this polarization level requires S/N = 500. Additionally, 
exposures are needed at 8 position angles in order to construct 
the Stokes I, Q, and U fluxes (eight to eliminate systematics). 
Consequently, the required exposure time in total for this 
detection level, not counting overhead, will be a little over 30 
minutes using a 1 -meter telescope. Although this exposure 
estimate is a lower limit (owing to neglect of background, 
inefficiencies, and telescope overheads), the required exposure 
is about 10% of the duration of the photosphere crossing, 
even smaller for the bulk of the circumstellar scattering enve- 
lope, and can be reasonably obtained with modest facilities 



equipped with polarimetric instrumentation. 

Although the original goal of microlensing surveys was 
to deduce the properties of dark matter in the Milky Way, it is 
clear that a vast range important of by-products have resulted 
from the survey effort, from catalogs of variable stars to obser- 
vations of finite source effects (as described in the Introduc- 
tion). Our contribution to the topic of finite source effects has 
been to point out how novel and valuable information about 
circumstellar envelopes might be obtained through polarimet- 
ric monitoring of events involving binary lenses and sources 
that may have substantial winds. 

Certainly, our models include some simpUfying assump- 
tions, such as ignoring polarisation from the photosphere 
and the effect of limb darkening. Neither of these are severe; 
indeed, both will tend to increase the peak polarisations above 
those predicted by our models. Photospheric polarisation is 
expected to be smaller than the circumstellar contribution for 
stars with significant winds, but its contribution should add 
to the Q and U fluxes constructively, and not destructively. 
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Limb darkening will tend to mollify the effects of the finite 
depolarisation factor, thereby increasing the peak polarisation 
from the inner wind where the density of scatterers is larger 
(although limb darkening will have little or no impact in the 
case of significant central cavities). We have also not allowed 
for optical depth effects, by way of multiple scattering effects 
and extinction of the photospheric emission. We intend to 
consider these effects in a separate paper. However, substantial 
scattering optical depths will be more important for quite 
dense winds, like those of red supergiants and asymptotic 
giant branch stars. For the more common red giant stars, the 
optically thin assumption will be a good assumption. 
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APPENDIX A: OBTAINING (3 FROM THE 
ASYMPTOTIC DECLINE OF THE STOKES Q-FLUX 

Consider a source that lies inside the caustic and is passing 
out of it. (The arguments that follow also hold true when the 
source first approaches the caustic.) After the photosphere of 
the star has completely passed out of the caustic, the polar- 
isation achieves a strong positive (by our convention) peak 
value. This occurs because only the scattering envelope lies 
interior to the caustic and so is subject to the selective ampli- 
fication, whereas the photosphere is amplified by an approxi- 
mately constant value. As the event evolves, more of the en- 
velope transits the fold caustic, and the polarised signal at a 
given moment is given primarily by the highest value of the 
polarised flux along the caustic itself. The scale of this polar- 
isation is set by the scattering optical depth of the envelope, 
but the slope of the polarised flux light curve is determined 
by the density distribution of scatterers. Here, we derive this 
relation between the lens position and the dependence of the 
polarisation on (3. 

As the star moves farther along from the caustic, the 
scattered light is accurately described by a point source. 
Equation llOt indicates that far from the star, Q{p,a) ~ 
Qo (Rh/p)'^^^ cos 2a. The polarised flux will be given by 



D2 




D2 



dx A{a 



dx ■ 



dyQ(p,a) 

- ^ {x + x,f~y^ 
X I ^ (x + XiY + y'^ 

1 (/3+l)/2 

(Al) 



(x + xi)'^ + y 

Factoring out (x xi) and making a suitable change of vari- 
able, the integration in y can be evaluated numerically for any 



particular value of fi. For /3 an integer, analytic integration for- 
mulae will apply. Because the result is not critical for our con- 
cerns, we simply denote the result of the y integral as Y{0), 
giving 



dx/^/x 
{x + xiY' 



(A2) 



Using a substitution with z — y^x/xi, the integral will have a 
similar form as the one for y, but with a different dependence 
on (3. We define the result of this integration to be X(P), lead- 
ing to 



{fJ)Y{p). (A3) 



Observationally, the lens position is linear with time t, 
and so log _Fq = — ^^^-^ log t+log to, where the other factors 
have been collected into the variable to- Of course this is a 
power-law, and its slope is directly related to the value of (3 
for the density distribution of scatterers in the envelope. For 
example, many of our models employ /3 = 2, in which case 
the asymptotic polarised flux varies as Fq oc t^/^. It should 
be pointed out that this limiting behavior will only be achieved 
in the tail of the polarised light curve, somewhat following the 
polarimetric peak. 
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